Rectification and Phase Locking For Particles on Two Dimensional Periodic 

Substrates 
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We show that a novel rectification phenomena is possible for overdamped particles interacting with 
a 2D periodic substrate and driven with a longitudinal DC drive and a circular AC drive. As a 
function of DC amplitude, the longitudinal velocity increases in a series of quantized steps with 
transverse rectification occurring near these transitions. We present a simple model that captures 
the quantization and rectification behaviors. 
PACS numbers: 05.45.Ac, 05.60. Cd, 05.45.Pq, 05.40.-a. 



Particles transported through periodic potentials un- 
der an AC drive exhibit a remarkable variety of nonlin- 
ear behavior. One of the most intensely studied phe- 
nomena found in these systems is the phase-locking that 
occurs when the external AC frequency matches the fre- 
quency of motion over the periodic potential, such as 
the well-known Shapiro steps observed in the V(I) curves 
of Josephson-junction arrays [||. Phase locking steps in 
transport curves are also seen in other systems such as 
charge density waves 1^ and vortices in superconductors 
with periodic substrates If the potential through 

which the particle moves is asymmetric, a rectification or 
ratchet effect can occur [|6[-p^ in both thermal |^,^ and 
deterministic ratchets H7 The ratchet effect has been 
studied in the context of biological motors and can 
be utilized for practical applications such as particle seg- 
regation atom transport in optical lattices 0, and 
fluxon removal in superconducting devices In de- 

terministic ratchets, as a function of AC amplitude, a 
complex series of phases appear in which rectification in 
either direction and current reversals can occur |^ . 

In this work we study a novel form of rectification and 
phase locking for an overdamped particle without ther- 
mal effects moving in a 2D periodic substrate under an 
applied longitudinal DC drive (f^c) and two AC drives: 
f%Q in the longitudinal direction 90 degrees out of phase 
from in the transverse direction. For AC amplitude 
A large enough that the stationary {ioc = 0) particle or- 
bit encircles at least one substrate potential maximum, 
we find that the longitudinal velocity increases in a se- 
ries of steps as toe increases, and that near the transition 
between two steps a rectification in the transverse direc- 
tion occurs. The steps correspond to drives at which the 
particle motion forms orbits commensurate with the sub- 
strate period. The transverse rectification occurs when 
the DC drive imposes an asymmetry on the particle orbit, 
interfering with the commensuration and allowing parti- 
cles to "leak" through the substrate in the transverse 
direction. The rectification is predominantly in one di- 
rection; however, reversals of the rectification also occur. 
We capture the qualitative features of the steps and rec- 
tification with a simple model, and show that our results 



are generic. The behavior we describe can be observed 
in vortices in superconductors with periodic pinning, as 
well as overdamped charged particles such as colloids in- 
teracting with periodic substrates. In addition, the rec- 
tification could be useful for particle species segregation, 
such as in polydisperse colloidal systems. 

As a model for vortices in superconductors or colloids 
in solution, we consider an overdamped particle moving 
in 2D interacting with a repulsive periodic substrate ac- 
cording to the equation of motion: = f, + toe + ^ac = 
rjVi, with rj = 1. The force from the substrate, a square 
array of side a, is f, = — V?7(r); the form of U{r) is 
discussed below. We consider a system of size 8a x 8a. 
The DC drive iuc is applied along the symmetry axis 
of the pinning array, in the x or longitudinal direction. 
The AC drive is ^ac — v4sin(Li;^t)x-|-i?cos(a;Bt)y. Note 
that there is no DC driving component in the y or trans- 
verse direction. We fix wa/wb = 1.0 and A — B, and 
examine both the longitudinal time averaged particle ve- 
locity < Vx > and the transverse velocity <Vy>. inc 
is increased from to 1.0 in increments of 0.0001, with 
3 X 10^ time steps spent at each drive to ensure a steady 
state. 

To model specific systems, we consider periodic sub- 
strate potentials created by pinned particles, such as vor- 
tices in a periodic array of holes or magnetic dots . 
Once all the holes are filled with a vortex, any remain- 
ing vortices sit in the interstitial regions between holes. 
An unpinned vortex experiences a smooth periodic sub- 
strate created by the interactions with the pinned vor- 
tices. Phase locking of AC and DC driven interstitial 
vortices can occur, and has been observed in experiments 
Q and simulations For {ac = 0, there is a finite 
DC depinning threshold for the interstitial vortices. For 
toe = 0, two crossed AC drives with small amplitudes 
A cause the vortex to move in a circle in the interstitial 
regions. At increasing A, there are stable vortex orbits 
which encircle one pinned vortex, then four, nine, and so 
on. We concentrate on A large enough to generate orbits 
encircling one or more pinned vortices. For the pinned 
particles, we use the potential for vortices in a thin film 
superconductor, U{r) = — ln(r), and employ a 
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FIG. 1. (a) Longitudinal velocity Vx vs driving force fnc 
which is applied in the x direction, (b) Corresponding trans- 
verse velocity Vy vs foe- 

summation technique [l^ for the long range interac- 
tion. We have also considered potentials for unscreened 
or screened charges with interaction of U{r) = 1/r and 
U{r) — e^'""/r, respectively. 

In Fig. 1(a) we show the longitudinal velocity Vx/{auj) 
versus {dc for a particle in a system with U{r) ~ ~ hi(7'), 
and A ^ B = 0.36. For these parameters, at ioc = 
the particle encircles one pin for 0.29 < A < 0.4. As 
shown, Vx increases in quantized steps with step heights 
of AVx — aoj. We label a step n according to the value of 
Vx on the step, Vx — nauj. The widths of the integer steps 
varies with A and lo in an oscillatory manner. There are 
also some fractional steps with heights AT4 — {pl<l)a'^ 
where p and q are integers such that p/ q < 1.0. 

In Fig. 1(b) we show the transverse velocity Vy/{aijj) 
versus foe- Since there is no net DC force in the trans- 
verse direction, Vy = Q iov most values of ioc- Near 
the steps in Vx, however, Vy is nonzero, indicating that 
rectification is occurring. The first rectification, in the 
positive y direction, occurs at the n = ^ 1 step, while 
at the n = 1 — > 2 step the rectification is in the negative 
y direction. For the higher steps, positive rectification 
regions with increasing widths appear, while near the 
n = 7 — > 8 step there are regions of negative rectifica- 
tion. 

Some of the rectification phases have very well defined 
heights of /SVy = Loa. These include the positive recti- 
fication phases at 3 ^ 4, 4 5, and 5 — > 6, as well as 
the negative rectification phases at 6 — > 7 and 7^8. 
On these rectification plateaus, we find that the particle 
moves in only one type of orbit. The other rectification 
phases do not have a well defined height, including 0—^1, 
2 — > 3, as well as some portions of the 6-^7 and 7^8 
regions. In these phases, for any fixed ^dc the particle 
jumps intermittently between different rectifying orbits 
with transverse velocities {p/q)auj. 

For all of the rectifying regions, if the polarity of the 
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FIG. 2. Trajectories for fixed foe for different 

non-rectyfing regions seen in Fig. 1. The black dots denote 
the location of the fixed particles or the potential maxima of 
the periodic substrate. Shown are steps with: (a) n = 0, (b) 
n = 1, (c) n = 2, (d) n — 4, (e) n — 5, and (f) n = 8. 

AC drive is reversed, Vx remains unchanged while Vy 
changes to —Vy. If the waiting time between DC drive 
increments is increased, the results are unchanged. The 
phases described here remain stable when the system is 
started from a fixed tjjc value, and are not transient 
phenomena. We have considered various system sizes 
with even and odd numbers of lattice constants a. As the 
system size varies, certain orbits become incommensurate 
with the system length and precess spatially; however, 
the velocity curves are not affected by the system size. 

In Fig. 2 we illustrate representative non-rectifying 
particle orbits for Vx steps of n = 0,1,2,4,5, and 8. 
For n — [Fig. 2(a)], Vx = 0, and the particle moves 
in a confined square orbit which encircles one pin. For 
n — 1 [Fig. 2(b)] there is a net motion in Vx, and the 
particle circles around one pin before moving over to the 
next plaquette. For n ~ 2 [Fig. 2(c)] the orbit does not 
encircle a pin but forms a small loop which repeats every 
second plaquette. For n — 4 and n — 5 [Fig. 2(d,e)], 
orbits similar to n = 2 occur, with the loop motion now 
repeating at every fourth or fifth plaquette, respectively. 
A similar process continues up through the n = 7 step. 
For the n — 8 step [Fig. 2(f)] and above, the particle 
is moving fast enough that the transverse width of the 
orbit is less than a, and no loops appear. 

Figure 3 shows representative rectifying orbits. The 
particle orbits differ below and above a given step, as il- 
lustrated in Fig. 3(a,b) for the n = 3 4 step. Below 
the step [Fig. 3(a)], the particle moves 3a in the x direc- 
tion and a in the positive y direction in a single period. 
A loop forms when the particle moves in the y direction. 
Above the step [Fig. 3(b)], the particle moves 4a in the x 
direction but still only a in the y direction in one period; 
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FIG. 3. Trajectories for fixed foe for different rectify- 
ing regions seen in Fig. 1. The black dots denote the loca- 
tion of the fixed particles or potential maxima of the peri- 
odic substrate, (a) right before n = 3 ^ 4 step, (b) right 
after n = 3 — > 4 step (c) right before n = 4 — > 5 step, 
(d) negative rectification phase below the n — 6 ^ 7 step 
(foe = 0.454), (e) negative rectification region on the n = 7 
step (foe = 0.465), (f) negative rectification region on n = 8 
step (foe = 0.514). 

therefore, Vy does not change at the 14 jump. In 
Fig. 3(c), below the n = 4 ^ 5 step, an orbit similar to 
that of Fig. 3(a) appears, except that the particle moves 
4a in the x direction in each period. The rectifying or- 
bit above the 4^5 step resembles that of Fig. 3(b), 
with the particle moving 5a in the x direction in one pe- 
riod. In Fig. 3(d) {foe = 0.454) we show the negative 
rectification phase below the n — 7 ^ 8 step. Here the 
particle jumps a in the negative y direction every seventh 
plaquette through a small kink. We do not observe loops 
in the trajectories for the negatively rectifying phases. 
Figure 3(e) {ioc = 0.465) illustrates trajectories for a 
particle in the negative rectification region close to the 
n — 7 ^ 8 step. The particle does not move in a spe- 
cific orbit, but jumps intermittently over time between 
different orbits with Vy = (j)/q)au)] however, Vx remains 
fixed. Similar intermittent trajectory patterns appear 
in the rectifying phases near n = 0^1,l— >2, and 
2 — > 3. Intermittent patterns occur only on the lower 
steps; above the n = 8 step, only stable rectifying trajec- 
tories occur, as illustrated in Fig. 3(f) {foe = 0.514). 

Rectifying phases occur for any A large enough that 
the particle trajectory at foe — encircles more than one 
pin. We have measured Vy/ aoj for the same system shown 
in Fig. 1, but with an orbit at ^4 = 0.42 that encircles four 
pins. As in Fig. 1, steps in Vx and transverse rectification 
in Vy appear, but the step heights are now A14 = 2aa;. 
We find that as A is further increased, orbits that stably 
encircle pins, with p integer, produce steps of height 
^Vx — paLu. We obtain very similar results for substrates 



with a l/r or e~'^'^/r interaction. 

The quantization of the particle velocity on the 
plateaus and the intermittent transitions between orbits 
can be explained using general properties of nonlinear 
maps. Define a map {x,y) — > {x' + nxa,y' + UyO), from 
the position of the particle at the start of a period to that 
at the end, where we may restrict to < x,y,x',y' < a, 
with nx,ny integer. If there is a stable fixed point, 
(x, y) — {x' , y'), then the particle translates by (ri^^a, nyo) 
in time uj~^ and so has average velocity Vx^ Vy quantized 
in multiples of aw, as found above. If the q — th power of 
the map has a stable fixed point, there are instead steps 
of fractional heights {p/q)auj. 

As {dc increases, the periodic orbit becomes unsta- 
ble, and a different periodic orbit with larger Vx appears. 
This new orbit will be the next stable periodic orbit at 
higher drive. The transition to the new orbit can occur 
in one of three ways. both periodic orbits are stable 

simultaneously, the particle velocity will depend on the 
initial conditions in the transition regime. This was not 
observed. (2) The second periodic orbit could become 
stable at the same time that the first orbit becomes un- 
stable. This behavior, which gives rise to infinitely sharp 
jumps in Vx, is not generic and hence not expected. (3) 
There can be a finite range of drive containing no sta- 
ble periodic orbits. Over this range, the average velocity 
is not quantized. If, however, some orbits are close to 
stable, the particle will spend long times in these orbits, 
giving rise to intermittent behavior. This behavior is 
consistent with what we observe. 

We now turn to a specific toy model illustrating some 
of these ideas. Consider a particle in a lattice of repulsive 
sites with a = 1, where the potential minima between re- 
pulsive sites are at integer x and y values. The y position 
of the particle is constrained to take only integer values, 
but the X position can be any real value. To model the 
translation of the particle through the lattice, we sepa- 
rate the X and y motion, so that the particle moves first 
(i.) right, then (ii.) down, then (iii.) left, then (iv.) 
up. (i): We apply the rule x ^ x + Vr- (ii.): If x is 
within 0.25 of an integer, x is set to that integer and y is 
decremented by one. (iii.): Apply x ^ x — vi. (iv.): As 
in (ii.) except y is incremented by one. Here Vr and vi 
are the velocity of the particle in the rightward and left- 
ward parts of the cycle, respectively. In steps (ii.) and 
(iv.), the particle will only move to a new y position if 
it reaches the minima between sites at the correct phase 
of the driving period, when transverse motion is possi- 
ble. In this case, the particle slips into the next row and 
the X coordinate of the particle is set to midway between 
the pinning sites. In Fig. ^ we show the time-averaged 
velocities Vx and Vy obtained with this model for fixed 
vi = 1.15 and increasing v^, representing increasing toe- 
As shown, this simple model produces both plateaus and 
ratcheting behavior. The sharp jumps in the velocity 
values are due to the discontinuity of the map function, 
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FIG. 4. (a) Time averaged longitudinal velocity Vx ob- 
tained from model, (b) Corresponding transverse velocity Vy. 

as in case (2) above; for smoother map functions, these 
jumps acquire a small but finite width. More compli- 
cated maps produce richer behavior, including occasional 
regimes of negative Vy. 

The ratcheting behavior in both the model and simu- 
lations occurs near transitions in Vx when the number of 
pinning centers the particle passes in one period changes, 
making it possible for the particle to interact asymmet- 
rically with the pinning sites. For a clockwise orbit, the 
particle moves rapidly on the upper portion of the orbit, 
and is likely to scatter off the pinning site below when the 
orbit does not quite match na. On the lower part of the 
orbit, however, the particle is moving more slowly, and 
is likely to slip between the pinning sites above in spite 
of a small mismatch. The particle thus tends to ratchet 
in the positive y direction. If the DC drive is reversed, 
downward motion should be preferred, as we observe. 

Finally, we note that much of this behavior is specific to 
two or more dimensions. Consider a map x x' , subject 
to X + a —>■ x' + a and dx'/dx > 0, true for overdamped 
motion in one dimension. It can then be shown that it is 
not possible to have stable periodic orbits with different 
values of the current. Systems with Shapiro steps do not 
exhibit jumps. 

The phases we have described should be experimen- 
tally observable for vortex motion in superconductors 
with periodic pinning arrays and Josephson-junction ar- 
rays for low flux filling or rational filling fractions where 
vortex-vortex interactions are reduced. The rectifying 
properties discussed here can also be a powerful tool for 
species separation. For example, colloids with different 
charges or DC mobilities driven through a periodic array 
would be active ratchets since the rectified velocities can 
be controlled by ncua, unlike thermal ratchets which rely 
on Brownian motion. Another system in which the longi- 
tudinal velocity steps and transverse rectification could 
be observed is electrons undergoing classical cyclotron 
orbits in anti-dot arrays |13| for orbits where electrons 



encircle at least one anti-dot. 

In conclusion, we have found that a novel form of rec- 
tification and phase locking occurs for a particle with a 
DC drive in the longitudinal direction and circular AC 
drive moving in a periodic substrate where the ampli- 
tude of the AC drive is large enough that the particle can 
encircle at least one substrate potential maximum. We 
find that the longitudinal velocity increases in a series of 
steps of height moa. Along these steps the particle moves 
in orbits that are commensurate with the periodicity of 
the substrate. For regions near the longitudinal step, 
transitions to rectification in either transverse direction 
occurs. We have specifically demonstrated this model 
for vortices in superconductors with periodic pinning ar- 
rays and overdamped charge particles. With a simple toy 
model we have shown that the qualitative features of the 
phase locking and rectification can be captured. 
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